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Abstract 

We present a new measurement of the time-dependent CP-violating parameters in B° — > ir + ir~ 
decays with 275 x 10 6 BB pairs collected with the Belle detector at the KEKB asymmetric- 
energy e + e~ collider operating at the T(45) resonance. We find 666 ± 43 B° — > ir + ir~ events 
and measure the CP-violating parameters: S n7T = —0.67 ± 0.16(stat) ± 0.06(syst) and A n7T = 
+0.56 ± 0.12(stat) ± 0.06(syst). We find evidence for large direct CP-violation with a significance 
greater than 4 standard deviations for any S wn value. Using isospin relations, we obtain 95.4% 
confidence intervals for the CKM angle <\> 2 of 0° < 4> 2 < 19° and 71° < 4> 2 < 180°. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 
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One of the unresolved mysteries in particle physics is the origin of CP violation. 
Kobayashi and Maskawa (KM) pointed out in 1973 that CP violation can be incorpo- 
rated as an irreducible complex phase in the weak-interaction quark mixing matrix in the 
standard model (SM) framework pj. The KM model predicts CP- violating asymmetries in 
the time-dependent rates of neutral B meson decays to the CP eigenstate 7r + 7r~ Q. In the 
decay chain of T(45) -> B°W -> (tt+tT ~)(/tag), one of the neutral B mesons decays into 
7r + 7r~ at time t nn and the other decays at time t tag to a final state / tag that distinguishes its 
flavor. The time-dependent decay rate is given by 

-|At[/r B o 

7^ (At) = — [1 + q ■ {S n « sin(Am d At) 

At b o 

+A mc ooa(Am d At)}], (1) 

where At = t nn — t tag , r^o is the B° lifetime, Am^ is the mass difference between the two 
neutral B mass eigenstates, and q — +1 (—1) for / tag = B° (B°). We measure S n7T and 
A n7T , which are the mixing- induced and direct CP- violating parameters, respectively. In the 
case where only a b — > u tree transition contributes to the decay -B — > n + 7i~ we would 
have S n7T = sin20 2 and A wn = 0. Because of possible contributions from b — > d penguin 
transitions that have different weak and strong phases, S nn may deviate from sin 202, and 
direct CP violation, A^ ^ 0, may occur. Our previous measurement based on a 140 fb _1 
data sample indicated large S n7T and A nn values |4j, while no significant CP asymmetry 
was observed by the BaBar Collaboration 0. It is therefore important to measure the 
CP- violating parameters with larger statistics. 

The measurement in this Letter is based on a 253 fb _1 data sample containing 275 x 
10 6 BB pairs collected with the Belle detector at the KEKB e + e~ asymmetric-energy (3.5 
on 8 GeV) collider [f| operating at the T(45) resonance. The T(4S) is produced with a 
Lorentz boost factor of (3^ = 0.425 along the z axis, which is antiparallel to the positron 
beam direction. Since the two B mesons are produced nearly at rest in the T(4S*) center-of- 
mass system (CMS), the decay time difference At is determined from the distance between 
the two B meson decay positions along the z-direction (Az): At = Az/cfl'j, where c is the 
velocity of light. 

The Belle detector |7] is a large- solid- angle magnetic spectrometer that consists of a 
silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to detect K\ mesons and to identify 
muons (KLM). A sample containing 152 x 10 6 BB pairs (Set I) was collected with a 2.0 cm 
radius beampipe and a 3-layer silicon vertex detector, while a sample with 123 x 10 6 BB 
pairs (Set II) was collected with a 1.5 cm radius beampipe, a 4- layer silicon detector, and a 
small-cell inner drift chamber 

We reconstruct B° — > tt + tt~ candidates using oppositely charged track pairs that are 
positively identified as pions by combining information from the ACC and the CDC dE/dx 
measurements. The pion detection efficiency is 90%, and 11% of kaons are misidentified 
as pions. We select B meson candidates using the energy difference AE = E* B — E^ eiim 



and the beam-energy constrained mass M bc = J (E^ cam ) 2 — (p* B ) 2 , where E£ eSLm is the CMS 
beam-energy, and E* B and p* B are the CMS energy and momentum of the B candidate. We 



define the signal region as 5.271 GeV/c 2 < M bc < 5.287 GeV/c 2 and \AE\ < 0.064 GeV, 
which corresponds to ±3 standard deviations (a) from the central values. 

We identify the flavor of the accompanying B meson from inclusive properties of particles 
that are not associated with the reconstructed B° — > 7r + 7r~ decay. We use q defined in Eq. ((*[)) 
and r to represent the tagging information. The parameter r is an event-by-event, Monte 
Carlo (MC) determined flavor-tagging dilution factor that ranges from r = for no flavor 
discrimination to r = 1 for unambiguous flavor assignment. It is used only to sort data 
into six r intervals. The wrong tag fractions for the six r intervals, wi (1 = 1,6), and the 
differences between B° and B° decays, Awi, are determined from data jlj Tol . 



To suppress the continuum background (e + e~ — > qq;q = u, d, s, c), we apply the technique 
used in Ref. 0. We form a likelihood ratio LR based on event topology variables and impose 
requirements on LR to suppress continuum events. The LR requirement is determined by 
optimizing the expected sensitivity using MC signal events and events in the sideband region 
in 5.20 GeV/c 2 < M bc < 5.26 GeV/c 2 or +0.1 GeV < AE < +0.5 GeV. As the tag-dilution 
variable r also suppresses continuum events, we optimize LR separately for each of the r 
bins. For all r bins we accept events having LR > 0.86; we then include additional events 
having lower LR thresholds of 0.50, 0.45, 0.45, 0.45, 0.45, and 0.20, respectively, for the six 
r bins. There are thus 12 distinct bins of LR-r for selected events. 

We extract 2,820 signal candidates by applying the above requirements and the vertex 
reconstruction algorithm used in Ref. [lOj to the data sample. Figure ^ shows the AE 
distributions for the events with (a) LR > 0.86 and (b) LR < 0.86 in the M bc signal 
region. The B° — > 7r + 7r - signal yield is determined from an unbinned two-dimensional 
maximum likelihood fit to the AE-M^ C distribution in the range of M bc > 5.20 GeV/c 2 
and —0.3 GeV < AE < +0.5 GeV with signal events plus contributions from misidentified 
B° — > K + tt~ events, the continuum background, and three-body B decays. We use a 
single Gaussian for the signal and B° —>■ K + tc~ events in AE and M bc . The continuum 
background shapes in AE and M bc are described by a first-order polynomial and an ARGUS 
function jllj], respectively. For the three-body B decay background shape, we employ a 
smoothed two-dimensional histogram obtained from a large MC sample. The fit to the 
subset with LR > 0.86 yields 415 ± 27 7r + n~ events and 154 ± 19 K + ir~ events in the signal 
region, where the errors are statistical only. The K + tt~ contamination is consistent with 
the K —>■ tt misidentification probability, which is measured independently. Extrapolating 
from the size of the continuum background in this fit, we expect 315 ± 3 continuum events 
in the signal region. We use MC-determined fractions as in |4J to calculate the numbers of 
decays for LR < 0.86. We expect 251 ± 16 tt+tt", 93 ± 12 K + n~ and 1, 592 ± 15 continuum 
events in the signal region. The contribution from three-body B decays is negligibly small 
in the signal region. 

We determine «S W7r and A nn by applying an unbinned maximum likelihood fit to the 
distribution of proper-time difference At. The probability density function (PDF) for 
the signal events is given in Eq. (JTJ) modified to incorporate the effect of incorrect fla- 
vor assignment w\ and Avj[. The distribution is convolved with the proper time inter- 
val resolution function R sig (At) in order to take into account the finite position resolu- 
tion [13 Q. The PDF for B° K+n- is V^{At, w h Awi) = (l/4r B o) e -l A *l/ r so [1 - 
qAwi + q(l - 2wi)Af n cos{Am d At)]. We use Af n = {A K n + A £ )/{1 + A K nA e ), where 
Artt = —0.109 ± 0.019 is the measured direct CP- violating parameter in B° — > K + tt~ 
decays ^3| , and A £ is the difference in the product of the pion efficiency and kaon misiden- 
tification probability between n + (K~) and ti~(K + ) divided by their sum 141 . The inclu- 
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FIG. 1: AE distributions in the M DC signal region for B 
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sion of A E changes the Ak-k value by 11%. We make use of the same resolution function 
i2 S ig(At) for the B° — > K + tt~ events. The PDF for the continuum background events is 
V q q(At) = 1/2 • (1 + qAqq) [(f T /2r qq )e^ A *l/ r ^ + (1 - f T )5(At)}, where f T is the fraction of the 
background with effective lifetime r qq , and 5 is the Dirac delta function. We use A qq = 
as a default. A fit to the sideband events yields A qq = +0.01 ± 0.01 (—0.00 ± 0.01) for 
the data in Set I (II). This uncertainty in the background asymmetry is included in the 
systematic error for the S n7T and A nn measurement. The background PDF V qq is convolved 
with a background resolution function R qq . All parameters in V qq and R qq are determined 
from sideband events. 

We define a likelihood value for each (ith) event as a function of S n7T and A^: 

/+oo 
[{fZV^t'^ h AwuS^A^) 
-oo 

+f%„V 9 K „(At', w h A Wl )}R^{AU - At') 
+f%V qq (At')R qq (At i - At')]dAt' 

+/oiPoi(Ati). (2) 

Here, the probability functions f™(k = tctt, Kit, or qq) are determined on an event-by-event 
basis as functions of AE and Mb c for each LR-r bin (m — 1, 12). A small number of signal 
and background events that have large values of At is accommodated by the outlier PDF, 
V i, with a fractional area / Q i. In the fit, S nn and A nn are the only free parameters and are 
determined by maximizing the likelihood function £ = Hi Pi, where the product is over all 
the B° — > 7r + 7r~ candidates. 

The unbinned maximum likelihood fit to the 2,820 B° — > 7r + 7r~ candidates containing 
666 ± 43 7T+7T- signal events (1,486 B° tags and 1,334 B° tags) yields = -0.67 ± 
0.16(stat) ± 0.06(syst) and A n7T = +0.56 ± 0.12(stat) ± 0.06(syst). The correlation between 
S n7T and A nw is +0.09. In this Letter, we quote the usual fit errors from the likelihood 
functions, called the MINOS errors, as statistical uncertainties 0]. Figures |2f a) andEfb) 
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(a) 470 candidates with q = +1, (b) 414 candidates with q = — 1. Raw asymmetry, Acp, in each 
At bin with (c) < r < 0.5 and (d) 0.5 < r < 1.0. The solid lines show the results of the unbinned 
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show the At distributions for the 470 B°- and 414 5°-tagged events in the subset of data with 
LR > 0.86. We define the raw asymmetry Acp in each At bin by Acp = (N + — N-)/ (N + + 
AL), where iV +( _) is the number of observed candidates with q = +1 (—1). FiguresE^c) and 
Etd) show the raw asymmetries for two regions of the flavor-tagging parameter r. 

The main contributions to the systematic error are due to the uncertainties in the vertex 
reconstruction (±0.04 for S n7T and iaoi f° r A^) and event fraction (±0.02 for S nn and 
±0.04 for Aw)', the latter includes the uncertainties in A q g and final state radiation. We 
include the effect of tag side interference on S nw (±0.01) and A^ 

(iaol)- Other sources 

of systematic error are the uncertainties in the wrong tag fraction (±0.01 for «S W7r and ±0.01 
for A wn ), physics parameters (r B o, Am d and Ak-k) (< 0.01 for and ±0.01 for A^), 
resolution function (±0.04 for S nn and ±0.01 for A n7T ), background At shape (< 0.01 for 
S n7T and < 0.01 for A n7T ), and fit bias (±0.01 for S wn and ±0.01 for A^). We add each 
contribution in quadrature to obtain the total systematic error. 

We carry out a number of checks to validate our results. The B° lifetime is measured with 
B° — > 7r + 7r~ candidates. The result is t b o = 1.50±0.07 ps, consistent with the world average 



value [17J. The CP fit to the sideband events yields no significant asymmetry. We check the 
measurement of Ann using a time-integrated fit and obtain A^ = ±0.52 ± 0.14, consistent 
with the time-dependent fit result. We also select B° — > K + n~ candidate events with 
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FIG. 3: AE distributions in the My, c signal region for the B° — ► tt + tt candidates with LR > 0. 
and 0.5 < r < 1.0 for (a) q = +1 and (b) q = — 1. 



charged tracks positively identified as kaons that have a topology similar to the B° — ► 7r + 7r~ 
signal events. The CP fit to the 4,293 B° — > K + ir~ candidates (2,207 signal events) yields 
Sk-k = +0.09 ± 0.08, consistent with zero, and Ak-k = —0.06 ± 0.06, in agreement with the 
world average value |lj|. With the K + ir~ sample, we determine t#o = 1.51 ± 0.04 ps and 
Airirf = 0.46 ± 0.03 ps -1 , which are also in agreement with the world average values [171 ] . 

To determine the statistical significance of our measurement, we apply the frequentist 
procedure described in Ref. that takes into account both statistical and systematic errors. 
The hypothesis of CP symmetry conservation, S n7r = A nw = 0, is ruled out at a confidence 
level (C.L.) of 1— C.L.= 5.6 x 10 -8 , equivalent to a 5.4(7 significance for one-dimensional 
Gaussian errors. The case of no direct CP violation, A nn = 0, is also ruled out with a 
significance greater than 4. Oct for any S nn value. 

Figure H3 shows the AE distributions for B° — > 7r + 7r~ candidates with LR > 0.86 and 
0.5 < r < 1.0 for (a) q = +1 and (b) q = — 1 subsets in the Mb c signal region. An unbinned 
two-dimensional maximum likelihood fit to the q = +1 (q = — 1) subset yields 107 ± 13 
(69 ± 11) 7T+7T-, 42 ± 9 (43 ± 9) K + n~ and 38 ± 1 (38 ± 1) continuum events in the signal 
box. The K + tc~ and continuum background yields are consistent between the two subsets 
as expected, while the tt + tt~ yields are significantly different; thus large direct CP-violation 
in B° — > 7r + 7r~ decays is manifest in the contrast of the two subsets. These results also 
support the expectation from SU(3) symmetry that A n7T ~ —3 ■ Ak-k (3- 

We constrain the ratio of the magnitude of the penguin to tree am plit udes \P/T\ and the 
strong phase difference 5 = 5p — 5t by adopting the notation of Ref. (l9j, where Sprr) is the 
strong phase of the penguin (tree) amplitude. By using (f>i = 23.5 ± 1.6° |l3T |. we find 95.4% 
confidence intervals of \P/T\ > 0.17 and -180° < 5 < -4°. 

To constrain (j) 2 , we employ isospin relations i2Q| and the approach of Ref. [ilj for the 
statistical treatment. We use the measured branching ratios of B° — > 7r + 7r~, 7r°7r° and 
B + — > 7T + 7r°, and the direct CP-asymmetry for B° — > ti it as well as our measured 
values of S n7T and A n7r taking into account their correlation. Figure EJ shows the obtained 
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FIG. 4: Confidence level as a function of the CKM angle 4>2 obtained with an isospin analysis using 
Belle measurements of S nn and A n7T . The dotted line indicates C.L. = 95.4%. 



C.L. as a function of 02- We find an allowed range for 2 a t 95.4% C.L. of 0° < <p2 < 19° 
and 71° < 2 < 180°. 

In summary, we have performed a new measurement of the CP- violating parameters in 
B° — > 7r + 7r~ decays using a 253 fb _1 data sample. We obtain S nn = —0.67 ± 0.16(stat) ± 
0.06(syst) and A W7r = +0.56 ± 0.12(stat) ± 0.06(syst). We rule out the CP-conserving case, 
S W7r = A W7r = 0, at the 5.4<r level. We find compelling evidence for direct CP asymmetry 
with 4. Oct significance. The results confirm the previous Belle measurement of the CP- 
violating parameters as well as the earlier evidence for direct CP violation in B° — > 7r + 7r~ 
decays |4j. 
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